WU-JIK-96-2 



Electrostriction of Polar Glasses 

J. I. Katz 

Department of Physics and McDonnell Center for the Space Sciences 
Washington U., St. Louis, Mo. 63130 
katz@wuphys. wustl. edu 

David R. Nelson 
Lyman Laboratory of Physics 
Harvard U., Cambridge, Mass. 02138 
(February 1, 2008) 

Abstract 

We develop a finite temperature theory for the susceptibility and electrostric- 
tion of isotropic substances in which permanent electric dipoles are restrained 
from free rotation by elastic forces. All parameters are constrained by the 
measured susceptibility and elastic constants. When applied to polyurethane, 
the predicted electrostriction is approximately consistent with some of the 
wide range of data. The saturation of the electrostriction at high field may be 
explained qualitatively if the dipoles consist of several amide groups locked 

together by crystallization of the hard segments of the polymer. 
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I. INTRODUCTION 



Zhenyi, et al. |l| discovered that the electrostriction of polyurethane film is large enough 
to be interesting. The measurements of Wang, et al. found a smaller electrostriction. The 
experiments are unexpectedly difficult, partly because the material is hygroscopic and its 
properties history-dependent, and partly because the measured electrostriction is affected 
by the confinement of the soft polyurethane by any stiffer material (such as electrodes) with 
which it is in contact. 

Comparatively little work exists ||[|] on the theory of electrostriction, and none directly 
applicable to this problem. We develop a simple theory of the electric susceptibility of 
polymers (such as polyurethane) with polar groups, and calculate electrostriction from the 
susceptibility. The theory has one parameter, which may be determined by comparing to 
the observed susceptibility. 

II. SUSCEPTIBILITY 

Polyurethane elastomers are block copolymers containing crystallizable and noncrystal- 
lizable segments M. The molecular structure of a generic polyurethane is shown in Figure 
[I], with x typically about 2. The crystallizable (hard) segments have attached amide groups 
with permanent electric dipole moments p ~ 3.6 x 10~ 18 esu-cm. Our model for the suscepti- 
bility considers the orientation of the dipoles in the electric field, and for the moment neglects 
the intrinsic susceptibility of the nonpolar portions of the polyurethane. For comparison, 
nonpolar polymers such as polyethylene typically have dielectric constants (resulting entirely 
from polarization of the chemical bonds) of about 2.3, while polyurethanes have dielectric 
constants of about 6.8. 

We make an independent particle approximation in which the dipoles in electrified 
polyurethane have two contributions to their energy — the usual electrostatic energy 

C/i = -p-E, (2.1) 
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and a term resulting from their rotation from their equilibrium orientations. This latter 
energy results from the elastic restoring forces of the polymer chains and of the bulk medium 
(entangled chains). For small deviations A9 <C 1 from the equilibrium position the elastic 
energy will, in general, be quadratic in AO. We therefore write 

U 2 = \s{A9)\ (2.2) 

where S is an elastic stiffness constant. For simplicity we assume that S is a scalar so that 
XJ 2 depends only on the magnitude of the angular deflection, and not on its direction; this 
may not be true for dipoles on a polymer, but in the absence of a detailed model for the 
orientation of the amide groups in polyurethane it is a reasonable approximation. If the 
applied electric field defines the 9 = axis of a polar coordinate system, and the orientation 
of the dipole defines the = plane and makes an angle 9 to the electric field, then A6 is 
obtained from 

cos AO = sin 9 sin 9 cos O + cos #0 cos 9, (2.3) 

where (6*0, <f>o) is the equilibrium position of the dipole in the absence of a field. 

Unfortunately 1273 becomes cumbersome when |2.3| is substituted, and the integrals re- 



quired to take thermal averages cannot be done analytically. Therefore, we replace |2T2 
by 

U 2 = S(l-cosA0). (2.4) 



To second order in AO 2.4 agrees with [2.2|. Because the actual calculation of U% for large 



A9 is a problem in nonlinear elasticity, which would be quite difficult, there is no reason to 



think that |2.4| is any less accurate than |2.2|. The form 2.4 may be rewritten 



U 2 = S - p • E, (2.5) 

where E is a vector with the dimensions of an electric field which points in the direction of 
the equilibrium position of the dipole (for E = 0), and which has magnitude E = S/p. The 



effects of the elastic forces on the dipole have now been incorporated into the field E, which 
will be assumed to be an isotropic random function of space with zero correlation length 



and fixed modulus E. The total energy U = U\ + U2 may be obtained from 2.1 and 2.5 



dropping the constant term in 2.5: 



£/ = -p-(E + E). (2.6) 

We ignore the electrostatic and elastic interaction between dipoles, in effect assuming 
their density to be low. The physics of randomly oriented dipoles embedded in an amorphous 
elastic polymer matrix is described by the random field E. If interactions between dipoles 
were to be included the problem would be that of a dipole glass, analogous to the familiar 
(but difficult) spin glass problem. 

We evaluate the mean polarization (p) in thermal equilibrium at temperature T. For a 
given dipole the energy |2]6| leads, by the usual elementary calculation, to a polarization of 
magnitude 

\(p) T \ = L(x c )p (2.7) 
along the direction of E + E, where the Langevin function 

L(x) = coth(x) - -, (2.8) 

x 

the combined normalized effective field strength is 

IE + El . , 

*c - LjA (2-9) 

and Qt denotes a thermal average only. 

The net polarization of the medium must, by symmetry, be in the direction (z) of the 
applied field E, and is found from 



T 



— J L(x c ) cosO'dtt, (2.10) 



p 47T 

where the overline represents a quenched angular average over all directions of E and 9' is 
the angle between E and E + E. Upon expanding L(x c ) and cos 9' in powers of S = E/E 



and £ = pE/kT = xoS, where xq = pE/kT, and keeping only terms of first order in E, we 
find 



^ = lsL(xo) + Ul'(x ). (2.11) 
p 3 3 

Note that while 5 1 and £ 1 are assumed, xq is generally not small. For a density n of 



independent dipoles the macroscopic polarization P = n(p)r may be obtained from [2.11J : 



E 

np— 
E 



2 1 
-L(x ) + -x L'(2;o) 



= npS/(o;o). (2.12) 
E 



The function /(x) is plotted in Figure |2|. The limits of [2.12j are 



-Ttor-, if x < 1; 
.fx„»l. 

The limit for xq <C 1 is a familiar elementary result, while the limit for x ^> 1 was given in 
slightly different form by Frolich ||. 

Polyurethanes typically have dielectric constants e w 6.8. Using the relation (in cgs 
units) e = 1 + 477%, this yields a susceptibility x — 0.46. The self-consistent molecular 
polarizability 7 (here a scalar, but in general a tensor) is 



P f <P*)t_^ 3 /e-l\ 

^ (jfcrj = ~E~ = 7 = 4^ UT2 J ' (2 ' 14) 

where the last equality is the Clausius-Mossotti equation || . This equation may be modified 
straightforwardly to allow for a dielectric constant e c of the continuum in which the dipoles 
are embedded. The result is 

(2 ' 15) 

For the polyurethane shown in Figure [I] with x = 2 and density 1.1 gm/cm 3 , the dipole 
density n = 3.4 x 10 21 cm -3 , assuming each amide group rotates independently. Then, with 
p = 3.6 x 10~ 18 esu-cm, the self-consistent solution to the transcendental equation |2.15| , 
taking T = 300°K and e c = 2.3, is pE/kT = 5.2, i.e. E = 6.1 x 10 4 cgs (= 1.8 x 10 7 V/cm). 
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It should be remembered that these numerical values depend on the assumed numerical 
values of n and p, which are poorly known. 

The twisting stiffness constant equivalent to the derived value of E is 5* = Ep = 2.2 x 
10~ 13 erg. This may be compared to the macroscopic elastic constants of polyurethane, 
which typically has a Young's modulus Y « 3 x 10 8 dyne/cm 2 . If the dipole is described 
by a rigid sphere of radius a, then rotation by an angle A9 implies a strain ~ AO, a 
tangential stress on the surface of the sphere ~ YA6, and a total torque ~ 2irYa 3 A6. The 
corresponding stiffness constant S ~ 2TrYa 3 , so that a ~ 5 A. This value is not unreasonable, 
although perhaps a bit large, but its significance is uncertain: The dipoles are mechanically 
coupled to their environment in complex ways. They are parts of a nearly inextensible 
and flexurally stiff covalently bonded polymer chain, which can shear quite easily with 
respect to its neighbors, except as limited by entanglements. These elastic properties are 
strongly temperature dependent, so that E will also be sensitive to temperature, at least 
for temperatures of order room temperature. Most of the temperature dependence of |2.11| 
and |2.12| is likely to come from the variation of E with T, rather than from the explicit 
dependence on T, except perhaps at very low temperatures. In fact, some polar polymers 
are known [|7| to have values of e which increase rapidly with T, as this argument would 
suggest. 



III. SATURATION 

The derivation of |2.12| is limited to small electric fields. The effective electric field E e ff 
acting on an object in an amorphous medium is (in cgs units) 

E eff = E + -P = e( — )^ 2.9E, (3.1) 



3 V 3 

where we have taken the empirical dielectric constant of polyurethane. In order for |2.12| to 



be valid the conditions pE e ff <C kT and E e ff -C E must both be met; when they are not 
met (p z ) becomes a significant fraction of p, the polarization saturates, and e, Xi an d the 
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electrostriction decrease below the values given by the low field theory. The available strong- 
field data QXJ] only describe the electrostriction, and show that it saturates at E « 10 7 V/m w 
300 cgs, corresponding to E e ff « 900 cgs, pE e ff/kT « 0.1, and E e ff/E ~ 0.02. We expect 
the susceptibility to saturate at the same field strengths as does the electrostriction, a 
prediction which is readily tested. 

The electrostriction saturates at fields substantially smaller than expected. This discrep- 
ancy may be resolved if the actual rotating dipole moments consist of a number of individual 
amide groups, locked together by crystallization of the hard segments of the polymer. Such 
a rigid block would have a value of p several times that of an individual amide group, and 
n would represent the density of such blocks and would be correspondingly smaller than 
the density of individual amide groups. The value of np would not be greatly changed by 
this locking of individual amides into a single giant dipole. Because f(x) is a slowly varying 
function for x ^> 1, E would also not be greatly changed. 



IV. ELECTROSTRICTION 

The Landau theory of electrostriction defines a free energy per unit volume of a body in 
an imposed stress field er^ and electric field Ef. 

F(aij, Ei) = fiu^ + ^Xu 2 mm - UijOij + ^(x^)tjPiPj + QijkiPiPjUM - EiPi, (4.1) 

where \i and A are Lame coefficients of elasticity, it^ is the strain tensor, is the stress 
tensor, {x~ l )ij is the inverse of the susceptibility tensor Xij defined by Pi = XijEj, Qijki 
is the electrostriction tensor, and Pi and Ei are the polarization and electric field. In this 
theory mj and P{ describe the internal state of the material, and characterize its response 
to the externally imposed parameters and E^. If we neglect fluctuations, which are small 



far from phase transitions, Uij and Pj may be found by minimizing 4.1 



By multiply differentiating the Legendre transform of 4.1 



G(<Tij ,Ei)=F + UijOi, + EiPi (4.2) 



in various orders we can relate Qijki to the derivatives of (x l )ij with respect to u^i- 

d d dG d d dG 
Qijkl ~ mdP^ ~ d^dPidR' (4 - 3) 



i.e., 



Qijkl = mw^ 1 = (44) 



where |C4] has been obtained by using the thermodynamic relations dG/duki = cr^i and 
(x~ r )ij — d 2 G / dPidPj. We then use our elementary microscopic model of x t° calculate 
Qijki from the second part of [CI] Q ijM = d(x~ 1 )ij/9uki- 



The application of a strain component M33 to an isotropic distribution of vectors E (the 
equilibrium orientations of the dipoles) rotates them. We assume that the magnitude E is 
not changed by the application of the strain field. We do not know if this is true, but it 
could not be described except by a completely phenomenological parameter, which we have 
no independent means of determining. The angle between E and the 9 = axis is rotated 
by an amount, found from elementary geometry 

A0 = -M33 sin 9 cos 6. (4.5) 

The calculation of susceptibility is then repeated, replacing 9 by 9 + A9, with the result, 
correct to first order in w 33 , 



733 = %^ = f[/(z )+ g(x )u 33 ], (4.6) 
&z E 

where the function g(x) is defined 

g ( x ) = -^[L(x)-xL'(x)}, (4.7) 
and is plotted in Figure 0. By considerations of symmetry (or explicit calculation) it is seen 



that jij = d(pi)T/9Ej is diagonal, as is Xiji so that (x ^33 = (X33) 1 - Then, to lowest order 
in M33, noting that An = —u 33 n and d{x~ 1 )ij I 'd^nj' 1 )^ = 1, we have 

n _ g(jT% _ E ( g{x )\ 

Q3333 ~ ~ ^tr { JW) r m 



We have neglected the dependence on strain of the susceptibility of the continuum medium 
in which the dipoles are embedded; any such dependence is likely to be very small, especially 
in a soft elastomer in which deformations are almost purely volume-conserving. 
The effects of strain components uu and U22 may be found similarly: 

A9 = (uu sin 2 (ft + u 2 2 cos 2 cf>) sin 6 cos 9, (4.9) 

where <fi is the azimuthal angle of E, measured from the y-axis. Then 

= *- = ^W) f 1 + w£) ■ (410) 

Similar calculations for U12, W13, and M23 (or symmetry arguments) show that 

Q 33kl = iork^l. (4.11) 

Similarly, because the distribution of E retains inversion symmetry in the strained state 
in all of these strained configurations Xij remains diagonal, so that Qijki = for i ^ j. 
Permutation of indices leads to the general result 



E 



Qijkl ' 



(l _ 9(xq) \ 
npf(x ) V f(x )J 



for i = j = k = I] 



E g(x ) \ „ . . . (4-12) 



1 H — — r for i = j 7^ k = I; 



npf(xo) V 2/(z )/ 
otherwise. 



We now calculate the electrostrictive strain produced by the application of an electric 

field and polarization in the z direction. Differentiating [4.1| with respect to Uki, and setting 
the derivative equal to zero in equilibrium yields 

= 2/MUki + S M Xu mm - a M + QijkiPiPj- (4.13) 

The electrostrictive term acts as an effective stress 

a' kl = -Q i3kl PiP 3 (4.14) 

= -QwuPl (4.15) 



Note that linearizing |4.15| about a state of nonzero E = (Eq + AE)z leads to 



Aoi 



ki 



-2Q 33kl ^^-P(E )AE; 



i.e., an effective piezoelectric response. Using |4.8| , |4.10| and [4. 1 1| yields 



npf(xo) 



2f(x )J 



(4.16) 



(4.17) 



E 



_ g(xo) \ p 2 . 

33 " npf(x )[ l f(x ))^ 



(4.18) 



the off-diagonal components are zero. Substitution in the elastic stress-strain relation (or 
use of |4.13l ) yields the strain 



Ull = «22 



+ 



E 



E gjxp) 
A/j,np f 2 (x ) ' 3Knpf(x ) 



Pi 



U33 



( E gjxp) | E 1_ 

\2fmpf 2 (x ) 3Knpf(x ) 



P 



3 ■ 



(4.19) 



(4.20) 



where the bulk modulus K = A + 



In order to rewrite |4.20| in terms of E 3 the Clausius-Mossotti equation |2.14j , p,15| may be 
solved for x m terms of f(xo), but because x is measured it is easier simply to substitute 
P = X E: 

1 



uu = ^ = ( E g(x ) + _E 

I Afmpf 2 (x ) 3Knpf(x ) i 



X E 3 



(4.21) 



U33 



/_^^(xo^ + _E 1_\ ^ 2 

\2finp f 2 (x ) 3Knpf(x )J 3 ' 



(4.22) 



For materials which are soft to shear, with \x K (as in many polymers) the second 



terms in parentheses of [4.20| and [4. 22| are small compared to the first terms. It is evident 
that materials with small shear moduli (or small Young's moduli, equal to 3/z in the limit 
fi/K — > 0) will show comparatively large electrostrictive strains. In addition, even though 
typically Q oc the strain u oc x an d materials with large x (but small Q) are usually 
the best electrostrictors IBfl. 
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V. COMPARISON TO EXPERIMENTS 



Substitution in 4.22 for PT6100S polyurethane 0, which has /x ~ 3 x 10 7 erg/cm 3 , 
taking K ^> /i, and using the estimated n and p and the value of E previously determined 
from the susceptibility, yields M33 = —7 x 1CT 9 cm 4 /esu 2 Ef , where the esu is the cgs unit 
of charge (1/(3 x 10 9 ) Coulomb). Zhenyi, et al. M find u 33 m — 2 x 10 -8 cm 4 /esu 2 E^, but 
with a range of values of a factor ~ 2 depending on the details of sample preparation and 
of mechanical constraint during the experiments. 

Wang et al. || find for 80CF-2 polyurethane M33 ~ — 3x 10~ 9 cm 4 /esu 2 E\ (approximately 
40% of the effect results from Maxwell stress rather than electrostriction). For this material 
/i « 1.0 X 10 s erg/cm 3 @ so that we predict (|4.22| ) m 33 « —2.4 x 10 -9 cm 4 /esu 2 Ef, 



approximately consistent with experiment. The result 4.22 contains no free parameters 



because the ratio E/(np) is determined from the measured dielectric constant by |2.15| , as is 
X, and \x is directly measured. 
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APPENDIX A: VOLUME TORQUES 

It is evident that there are no body forces on a uniformly electrified homogeneous mate- 
rial, but it is still of interest to examine the strain field u produced by a point torque in an 
elastic medium, such as the permanent dipoles in electric fields discussed in this paper, in 
order to see if this strain contributes to the electrostriction. This result not readily found 
in the elasticity literature. 
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Begin with Kelvin's result for the displacement field u produced by a point force F at 

the origin of an elastic medium ||: 

u(r) = ( IpL ) A3-4T)F + fl ( fl .F)A 
V ; V8ttF(1 -a) )\ r J' y 1 

where a is Poisson's ratio, r denotes the point at which the displacement is measured, and 
n is the unit vector in the direction of r. By adding to [Al| the displacement produced 
by a force — F acting at a point Ar we find the displacment produced by a point torque 
N = -Ar x F: 

/ N ( 1 + a \ N X r /AON 

"W = Ur(i- g ) J— (A2) 

A dipole p in an electric field E is subject to a torque 

N = p x E. (A3) 

In a medium with an inversion-symmetric distribution of dipole orientations (p) = and 
(N) = 0. Each dipole is rotated by an angle A.6 = ~N/S. 

The rotation of the dipoles changes the direction of their moments, and therefore the 
torque. The change in p is 

Ap^A^xp = ^P^ ' pX ^ Xp = Ep2 -P' p - E ) , (A4) 

and yields an incremental torque, second order in E (like the electrostriction): 

N = ApxE = - (pXE >< p - E > . (A5) 



Although [A5| is second order in both p and E, in any volume element it will average to 
zero for an inversion-symmetric distribution of p, and hence it is not a useful model of 
electrostriction. 



APPENDIX B: CORRELATION OF ELECTROSTRICTION WITH THERMAL 

EXPANSION 

Uchino, et al. ||ItJ note that the electrostriction coefficient Q and the thermal expansion 
coefficient a are empirically related by Q oc a 2 for a wide variety of materials, including one 
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polymer (PVDF), pyrex glass, ceramic piezoelectrics, and ionic crystals (see their Figure 1); 
there is only a factor of 2 (in a) or 4 (in Q) scatter about the power law fit, which extends 
over a range of about 30 in a (1000 in Q). This empirical relation has a very elementary 
explanation. 

The simplest model of thermal expansion is a particle in a weakly anharmonic potential 
V = ax 2 + bx 3 . If the particle has an energy kT, its turning points are shifted by Ax ~ 
—bkT/(2a 2 ), giving a linear thermal expansion coefficient a = \b\/(2a 2 a ), where a is the 
distance between the potential minimum and the origin (or the nearest neighboring particle). 

The analogous simple model of electrostriction considers two particles in the well, with 
charges ±q, and adds an electric field E. The particles are now in potentials V± = ax\+bx\±. 
qEx±, and their equilibrium (T = 0) positions are given by x± = —a + \J a? =F 3bqE/ (36). 
The displacement of their barycenter is Ax = (x + + x_)/2 ps —3bq 2 E 2 /(8a 4 ), giving an 
electrostrictive coefficient Q = — 36g 2 /(8a 4 a ) (note that b < for realistic potentials). 

If the chief variation between materials is in a then Q oc a 2 , as observed; the scatter 
about this relation arises from variations in b, a , and the effective charge q. The parameter 
a is related to the bulk modulus, which varies a great deal between materials; it is not 
surprising that variations in the other variables are less important, though not insignificant. 

APPENDIX C: ELECTROSTRICTION IN THE MELT 

It is possible to construct a very simple model of electrostriction of a melted poly- 
urethane, with no elastic resistance to rotation of the permanent dipoles. Suppose the 
polymer configuration to be a random walk which deviates from isotropy because an applied 
electric field distorts the equilibrium distribution of the orientations of the dipoles (for small 
anisotropy the self- avoiding property of the walk will not change the degree of anisotropy). 
Consider dipole components p± perpendicular to a fraction / of the units of the polymer 
chain, but otherwise free to rotate. The distribution of dipole orientations making an angle 
9 with respect to an electric field in the z direction is given by the elementary result 
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where v = p±E e ff/kT and E e ff = (1 + -fx)E * s the effective field at a dipole in a 
random medium. Each polymer r(s) can now be regarded as a biased random walk, 
whose mean square end-to-end distance grows anisotropically as a function of arc length 
([rj(s) — ri(0)][rj(s) — r,-(0)]) = (2DijS) 2ri . For chain elements of unit length the random 
walk is described by a diagonal diffusion tensor, correct to second order in u, 

D xx = \ ((y 2 d ) + (zl)) = I (l + (C2) 

Dyy = \ ((4) + (4)) = I (l + ^) (C3) 

^ = ^ (<*2> + <i/2>) = 5 (i - ^) , (C4) 

where the subscript d refers to the components of the dipole vectors. The corresponding 
strain is 

«--ii(i) 2 ( 1 + ^) 2 ^ — ^ «*> 

where we have taken / = 0.2, rj = 0.6, p± = 2.5 x 10~ 18 esu-cm, T = 300°K, and x = 0.46 
(corresponding to solid polyurethane with e = 6.8). This strain will make the material 
optically active, but in order to calculate this it would be necessary to know the intrinsic di- 
electric anisotropy of the polymer chain at optical frequencies; optical activity is determined 
by the strain alone, and is not specific to its electrostrictive origin. 
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FIGURES 

FIG. 1. Molecular structure of polyurethane elastomer. 
FIG. 2. The functions f(x) and g(x) 
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